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Abstract 

We present two integrable spin ladder models which possess a general free pa- 
rameter besides the rung coupling J. The models are exactly solvable by means of 
the Bethe ansatz method and we present the Bethe ansatz equations. We analyse 
the elementary excitations of the models which reveal the existence of a gap for 
both models that depends on the free parameter. 
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Spin ladder systems continue to attract attention motivated by experimental realisa- 
tions in quasi-one- dimensional systems [0. These materials display novel features and 
with the continued development of new systems, there has been an impressive amount of 
progress in the theoretical understanding of such systems. However, a greater flexibility 
through the introduction of tunable free parameters within the well established mathe- 
matical frameworks would be of considerable advantage and forms the main aim of the 
present work. 

It has been shown that ladder systems are reasonably well approximated by Heisenberg 
interactions, which involve bilinear exchanges || . While these models are not exactly solv- 
able, several more general systems have been proposed in which solvability is guaranteed 
through the use of an extension of the symmetry algebra [0-0] . There has also been the 
introduction of systems involving interactions beyond nearest neighbour exchanges which 
demonstrate remarkably interesting behaviour and also prove to be exactly solvable. For 
example, the general 2- leg spin ladder system with biquadratic interactions The 
physical importance of these types of interactions has been addressed in ||. 

Subsequently other generalised integrable spin ladders have been proposed 



As is well known, integrability facilitates the use of long established techniques in or- 
der to determine the physical properties of such systems. However, in these cases, no 
free parameters other than the rung coupling are present due to the strict conditions of 
integrability. 

In a recent article [15| an integrable model containing an additional free parameter was 



presented as a generalisation of the model presented by Wang [16]. In this instance, the 
free parameter plays the role of an anisotropy variable and it was shown that the critical 
value of the rung coupling which defines a Pokrovsky-Talapov phase transition between 



a gapped and gapless phase was dependent on this anisotropy. In [17], we note the study 



of a family of spin ladder Hamiltonians which also have free parameters, although in this 
case the construction has a different mathematical origin. It is clear that this is a topic 
that warrants further investigation, since the availability of tunable parameters yields a 
richer phase structure. 

In this article, we present two new integrable generalised spin ladders, based on the 
SU(1\3) and SU(3\1) symmetries, containing an extra parameter. The free parameter 



arises in the models as a special choice of the multiparametric versions ||18|| . The models 
are integrable in the sense that they contain an infinite number of conservation laws 
and can be derived from a solution of the Yang-Baxter equation. This property is also 
of physical importance as it provides a means to improve our understanding of general 
correlated systems (see for example [|19[). We present the Bethe ansatz solution from 



which the physical properties of the systems may be obtained. 

An important characteristic of ladder systems, both from a theoretical and experimen- 
tal point of view, is the quantum phase transition between gapped and gapless phases. 
The spin gap is vital for superconductivity to occur under doping, whilst from a mathe- 
matical perspective, the size of the gap is dependent on the relative strength of the rung 
interaction parameter. We address this issue as we analyse the ground state and first 
excited states of the models. Interestingly, we are able to show that for both systems a 
gap persists in the spectrum of the elementary excitations and indeed the gap depends 
on the extra parameter. 
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We begin by introducing the first generalised spin ladder model, for which the explicit 
global Hamiltonian is of the form 



N 



1 



3=1 

and the local Hamiltonians are given by 

hj+t = 7(1 + ^K+iX 1 + + W^ 7 +i + ^ +1 ){rfrf +1 + rr T / +1 

1 
2 



1 

+ a^^r 1 r/rr +1 + f r-r^ + V 1 + t aja+ +1 ){l + r/r/ +1 ) 



(1 + <)(i + ^ +1 )(i + t-/)(i + t-; +1 ). 



The parameters <fj and fj represent Pauli matrices acting on site j of the upper and lower 
legs respectively, J is the strength of the rung coupling that can take arbitrary real values 
and t is a general independent parameter. The number of rungs is denoted by N and 
periodic boundary conditions are assumed. 

The integrability of this model is assured by the Quantum Inverse Scattering Method 
20fl and by the fact that it can be mapped to the Hamiltonian given in equation (|2]) 



below. This Hamiltonian can be derived from an R— matrix obeying the Yang-Baxter 
algebra |2l| for J = 0, while for J ^ 0, the rung interactions take the form of a chemical 
potential term. We find that 
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where 
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a=0 
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+t 
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2XfXf +l . 



In the above, X^ = \aj) (/3j\ are the Hubbard operators with \aj) being the orthogonalised 
eigenstates of the local operator ((fj.fj). 

The i?-matrix we use is a special case of a more general multiparametric version. 
(For a general construction of multiparametric models, see For the purposes of the 
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present work, it is necessary to only retain one parameter. The i?-matrix is as follows, 
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and obeys the Yang-Baxter algebra [21] 



R\2{x - y)R 13 (x)R 2 3(y) = R 2 3(y) Ri 3 (x) R 12 (x - y). 



(3) 



(4) 



From this solution originates the Hamiltonian (||D for J = by the standard procedure |20[ 



P—R{x)\ x=0 , 
ax 



where P is the permutation operator. 

The model is exactly solvable by the Bethe ansatz method 
ansatz equations (BAE) are given by the expressions, 



and the resulting Bethe 
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The corresponding energy eigenvalues of the Hamiltonian (0) are, 
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+ 2J -(1 + 2J)JV, 



(6) 
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where Xj are solutions of the BAE 

From the Bethe ansatz solution, we can determine the behaviour of the ground state 
and elementary excitations of the system. The reference state becomes the ground state 
when the relation J > — 1 + ~(t + is satisfied. For N sites, the ground state energy is 
E = — (1 + 2 J) N, which in terms of the Bethe ansatz calculations, corresponds to the 
reference state characterised by Mi = M 2 = M 3 = 0. 

To describe an elementary excitation, we chose M 2 = M 3 = and Mi = 1 in the BAE 
which, from equation (B|), yields an energy expression of the form 

Ei = xr^nji + 2J - C 1 + 2 W ( 7 ) 

where A = | (fzi)- ^ * s apparent that there is a gap of 

A = 2(J + i-i(f + r 1 )). (8) 

In the limit t = 1, this solution corresponds to A — > oo indicating that a gap of A = 2J 
persists. We note that this is agrees with the suggested numerical and experimental results 
of spin ladder systems |]]. 

We move on to introduce the second integrable spin ladder model which also contains 
a free parameter. The global Hamiltonian reads 
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(2) 
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where 
k 
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o ^)(l-^ +1 )(l-r/)(l-r^). 

The exact solvability of the above Hamiltonian, as for the previous case, lies in the 
fact that it too can be mapped to a Hamiltonian given below by equation (10). Once 
again this Hamiltonian is derived from an i?-matrix solution of the Yang-Baxter algebra 
for J = 0, while for J ^ the rung interactions take the form of a chemical potential 
term. The Hamiltonian has the form, 
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For J = 0, the model is derived, in a similar manner as for the above case, from 
a multiparametric i?-matrix for which only one parameter is relevant for the present 
discussion. The i?-matrix is given by 
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and satisfies the Yang-Baxter algebra @ 
can be solved and the resulting BAE are, 



(11) 



Utilising the Bethe ansatz method this model 
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where Xj are solutions of BAE ([12]). 

For sites, the ground state is given by a product of rung singlets when 
J > 1 + |(t + t^ 1 ) an d the energy is So = (1 — 2J)N . This is in fact the reference state 
used in the Bethe ansatz calculations and corresponds to the case M\ = M2 = M3 = 
of the BAE (0). To describe an elementary spin-1 excitation, we choose M\ = 1 and 
M2 = M3 = in the BAE which gives the minimal excited state energy, 

£l = " AfTTTi + 2 J + (1 " 2 J) N ' (14) 
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where Ai = | (fqrj")- The energy gap can easily be calculated and is found to be 

1 



A 1 = 2^J-l-^ + r 1 )J. (15) 

The value J c = 1 + |(t + t -1 ) indicates the critical line at which the transition from 
dimerized phase to the gapless phase occurs. 

In conclusion, we have presented two new spin ladder models derived as special cases of 
multiparametric versions of 577(3|1) and S77(l|3) invariant solutions of the Yang-Baxter 
equation, maintaining one free parameter besides the rung coupling J. Upon investigation 
of the solutions of the BAE's to determine ground state and elementary excitations, we 
have shown that both models exhibit a gap that depends on the extra parameter. Our 
results show similar generic properties to the SU(4) model studied in j!5| and is very 



suggestive that such multiparametic extensions will, in general, always have an influence 
on the physical characteristics of these models, and in particular the critical value of the 
rung coupling. 
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